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• Abstract From the beginning of the subject, calculations of quantum vacuum ener- 

Z^ [ gies or Casimir energies have been plagued with two types of divergences: The total 

energy, which may be thought of as some sort of regularization of the zero-point en- 
ergy, ^ jh(0, seems manifestly divergent. And local energy densities, obtained from 
the vacuum expectation value of the energy-momentum tensor, (roo), typically di- 
_. , verge near boundaries. These two types of divergences have little to do with each 

flj ' other The energy of interaction between distinct rigid bodies of whatever type is 

(~| I finite, corresponding to observable forces and torques between the bodies, which 

can be unambiguously calculated. The divergent local energy densities near sur- 
faces do not change when the relative position of the rigid bodies is altered. The 
self-energy of a body is less well-defined, and suffers divergences which may or 
may not be removable. Some examples where a unique total self-stress may be eval- 

CO i uated include the perfectly conducting spherical shell first considered by Boyer, a 

perfectly conducting cylindrical shell, and dilute dielectric balls and cylinders. In 
these cases the finite part is unique, yet there are divergent contributions which may 

•Q I be subsumed in some sort of renormalization of physical parameters. The finiteness 

of self-energies is separate from the issue of the physical observability of the effect. 
The divergences that occur in the local energy-momentum tensor near surfaces are 
distinct from the divergences in the total energy, which are often associated with 
energy located exactly on the surfaces. However, the local energy-momentum ten- 

^^ ' sor couples to gravity, so what is the significance of infinite quantities here? For the 

H I classic situation of parallel plates there are indications that the divergences in the 

local energy density are consistent with divergences in Einstein's equations; corre- 
spondingly, it has been shown that divergences in the total Casimir energy serve to 
precisely renormalize the masses of the plates, in accordance with the equivalence 
principle. This should be a general property, but has not yet been established, for 
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example, for the Boyer sphere. It is known that such local divergences can have no 
effect on macroscopic causality. 



1 Introduction 

For more than 60 years it has been appreciated that quantum fluctuations can give 
rise to macroscopic forces between bodies |lj. These can be thought of as the sum, 
in general nonlinear, of the van der Waals forces between the constituents of the bod- 
ies, which, in the 1930s had been shown by London 15] to arise from dipole-dipole 
interactions in the nonretarded regime, and in 1947 to arise from the same inter- 
actions in the retarded regime, giving rise to so-called Casimir- Polder forces ||3]- 
Bohr m apparently provided the incentive to Casimir to rederive the macroscopic 
force between a molecule and a surface, and then derive the force between two con- 
ducting surfaces, directly in terms of zero-point fluctuations of the electromagnetic 
fields in which the bodies are immersed. But these two points of view — action at a 
distance and local action — are essentially equivalent, and one implies the other, not 
withstanding some objections to the latter ||5l. 

The quantum-vacuum-fluctuation force between two parallel surfaces — ^be they 
conductors or dielectrics IS 111 El — was the first situation considered, and still the 
only one accessible experimentally. (For a current review of the experimental situa- 
tion, see ||9][T0|) Actually, most experiments measure the force between a spherical 
surface and a plane, but the surfaces are so close together that the force may be 
obtained from the parallel plate case by a geometrical transformation, the so-called 
proximity force approximation (PFA) ifTTlfTSlfTSl . However, it is not possible to find 
an extension to the PFA beyond the first approximation of the separation distance 
being smaller than all other scales in the problem. In the last few years, advances 
in technique have allowed quasi-analytical and numerical calculations to be carried 
out between bodies of essentially any shape, at least at medium to large separation, 
so the limitations of the PFA may be largely transcended. (For the current status of 
these developments, see the contributions to this volume by Emig, Jaffe, and Rahi, 
and by Johnson; for earlier references, see, for example lfT4ll .) These advances have 
shifted calculational attention away from what used to be the central challenge in 
Casimir theory, how to define and calculate Casimir energies and self-stresses of 
single bodies. 

There are, of course, sound reasons for this. Forces between distinct bodies are 
necessarily physically finite, and can, and have, been observed by experiment. Self- 
energies or self-stresses typically involve divergent quantities which are difficult to 
remove, and have obscure physical meaning. For example, the self-stress on a per- 
fectly conducting spherical shell of negligible thickness was calculated by Boyer in 
1968 ifTSJI . who found a repulsive self-stress that has subsequently been confirmed 
by a variety of techniques. Yet it remains unclear what physical significance this 
energy has. If the sphere is bisected and the two halves pulled apart, there will 
be an attraction (due to the closest parts of the hemispheres) not a repulsion. The 
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same remarks, although exacerbated, apply to the self-stress on a rectangular box 
lfT6l [TTl [TSl [191 . The situation in that case is worse because (1) the sharp corners 
give rise to additional divergences not present in the case of a smooth boundary (it 
has been proven that the self-energy of a smooth closed infinitesimally thin conduct- 
ing surface is finite 1201 |2T1 ). and (2) the exterior contributions cannot be computed 
because the vector Helmholtz equation cannot be separated. But calculational chal- 
lenges aside, the physical significance of self-energy remains elusive. 

The exception to this objection is provided by gravity. Gravity couples to the local 
energy-momentum or stress tensor, and, in the leading quantum approximation, it 
is the vacuum expectation value of the stress tensor that provides the source term 
in Einstein's equations. Self energies should therefore in principle be observable. 
This is largely uncharted territory, except in the instance of the classic situation of 
parallel plates. There, after a bit of initial confusion, it has now been established 
that the divergent self-energies of each plate in a two-plate apparatus, as well as the 
mutual Casimir energy due to both plates, gravitates according to the equivalence 
principle, so that indeed it is consistent to absorb the divergent self-energies of each 
plate into the gravitational and inertial mass of each ll22l I23JI . This should be a 
universal feature. 

In this paper, for pedagogical reasons, we will concentrate attention on the 
Casimir effect due to massless scalar field fluctuations, where the potentials are de- 
scribed by 5-function potentials, so-called semitransparent boundaries. In the limit 
as the coupling to these potentials becomes infinitely strong, this imposes Dirichlet 
boundary conditions. At least in some cases, Neumann boundary conditions can be 
achieved by the strong coupling limit of the derivative of 5-function potentials. So 
we can, for planes, spheres, and circular cylinders, recover in this way the results 
for electromagnetic field fluctuations imposed by perfectly conducting boundaries. 
Since the mutual interaction between distinct semitransparent bodies have been de- 
scribed in detail elsewhere Il24l|25ll26l , we will, as implied above, concentrate on 
the self-interaction issues. 

A summary of what is known for spheres and circular cylinders is given in Table 
ffl 



2 Casimir Effect Between Parallel Plates: A 5-Potential 
Derivation 

In this section, we will rederive the classic Casimir result for the force between 
parallel conducting plates U. Since the usual Green's function derivation may be 
found in monographs ll38l . and was for example reviewed in connection with current 
controversies over finiteness of Casimir energies li36l . we will here present a differ- 
ent approach, based on 5-function potentials, which in the limit of strong coupling 
reduce to the appropriate Dirichlet or Robin boundary conditions of a perfectly con- 
ducting surface, as appropriate to TE and TM modes, respectively. Such potentials 
were first considered by the Leipzig group ||39ll40l . but more recently have been 
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Table 1 Casimir energy (£) for a sphere and Casimir energy per unit length (#) for a cylinder, both 
of radius a. Here the different boundary conditions are perfectly conducting for electromagnetic 
fields (EM), Dirichlet for scalar fields (D), dilute dielectric for electromagnetic fields [coefficient of 
(e — 1)^], dilute dielectric for electromagnetic fields with media having the same speed of light (co- 
efficient of (^^ = [(e — l)/(e + 1)]^), perfectly conducting surface with eccentricity 8e (coefficient 
of 8e^), and weak coupling for scalar field with 5-function boundary given by j60t . (coefficient of 
A^/a-). The references given are, to the author's knowledge, the first paper in which the results in 
the various cases were found. 

Type ^Sphere" ^^Cylinder"^ References 



EM 


+0.04618 


-0.01356 


nudii 


D 


+0.002817 


+0.0006148 


mm 


(e-l)2 


+0.004767 = 


_ 23 r, 
1536jr ^ 


mm 


^' 


+0.04974 = 


m 


02(331 


5e^ 


±0.0009 





Elllll 


X^/a' 


+0.009947 -- 


-^. 


mm 



the focus of the program of the MIT group BTl l42l l43l 1441 . The discussion here is 
based on a paper by the author 1451 . (See also B6l .) (A multiple scattering approach 
to this problem has also been given in IZSl .) 

We consider a massive scalar field (mass n) interacting with two 5-function po- 
tentials, one at X = and one at x = a, which has an interaction Lagrange density 

i^int = -\x5{x)<^\x) - ]^X'5{x-a)<^^{x), (1) 

where the positive coupling constants A and A' have dimensions of mass. In the 
limit as both couplings become infinite, these potentials enforce Dirichlet boundary 
conditions at the two points: 

A,A'^-: (^(O),0(fl)^O. (2) 

The Casimir energy for this situation may be computed in terms of the Green's 
function G, 

G{x,x')=i{T<^{x)<^{x')), (3) 

which has a time Fourier transform, 

G{x,x')^ /— e-"»('-'')^(x,x';a)). (4) 

J 2k 

Actually, this is a somewhat symbolic expression, for the Feynman Green's func- 
tion (O implies that the frequency contour of integration here must pass below the 
singularities in (O on the negative real axis, and above those on the positive real axis 
P7ll48l . Because we have translational invariance in the two directions parallel to 
the plates, we have a Fourier transform in those directions as well; 
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^(x,y ; co) = l ^e"'-('-"-')^^(x,x'; k), 



(5) 



where K^ ^ fi^ + k^- ap- 



The reduced Green's function in (|5]l in turn satisfies 

Kx,x')-5(x-x'). 



-4^ + ic2 + A5(x)+A'5(x-fl) 



(6) 



This equation is easily solved, with the result 



KV) = ^.-^I-I 



1 



A' 






2coshK'|x — x'l 



2k- V 2k-/ 2k- V 2k- 



(7a) 



for both fields inside, < x,x' < a, while if both field points are outside, a < x,x'. 



^M = ^e-l- 



1 



,-K-(x+V-2a) 



2K-Zi 

2k- V 2k-/ 2k- 



2k- 



^2Kcr 



(7b) 



For ji:,ji:' < 0, 



?{x,x') 



1 



,K-(x+y) 



2k- 



2K-Zi 

A 



2k- V 2k-/ 2k- 



2k- 



,2kci 



Here, the denominator is 



' 2K-y V 2k J (2k-)2 



(7c) 



(8) 



Note that in the strong coupling limit we recover the familiar results, for example, 

inside 

sinh K-x< sinh K-(x> — a) 



A, A' 



?ix,x') - 



(9) 



K sinh Ka 

Here x>, x< denote the greater, lesser, of x,x'. Evidently, this Green's function van- 
ishes at X = and at x = a. 

Let us henceforward consider jX ~ 0, since otherwise there are no long-range 
forces. (There is no nonrelativistic Casimir effect.) We can now calculate the force 
on one of the 5-function plates by calculating the discontinuity of the stress tensor, 
obtained from the Green's function ^ by 
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{Tn= (d''d"~\8^'d'd',y^Gix,x') 



Writing a reduced stress tensor by 



d« [ (dk) 
2nJ [Inf 



we find inside, just to the left of the plate at x = a. 



{-K^ + d,d,,)g{x,x') 



'2it^ ^2^Z 



(10) 

(11) 

(12a) 
(12b) 



From this we must subtract the stress just to the right of the plate at x = a, obtained 
from dTbb . which turns out to be in the massless limit 



f« 



K 

'2i' 



(13) 



which just cancels the 1 in braces in (Il2bb . Thus the pressure on the plate dXx = a 
due to the quantum fluctuations in the scalar field is given by the simple, finite 
expression 



^=(7:,v)|^^„_-(7:, 



1 



^="+ 32%- a^ 



dy/ 



1 



(y/(Afl) + l)(y/(A'«) + l)e.v-l' 

(14) 
which coincides with the result given in Il44l|49l . The leading behavior for small 
A = A' is 

^™^-^T^T' A<1, (15a) 

while for large X it approaches half of Casimir's result |[T] for perfectly conducting 
parallel plates, 



dTE 



A>1. 



(15b) 



480a4 ' 

We can also compute the energy density. Integrating the energy density over all 
space should give rise to the total energy. Indeed, the above result may be easily 
derived from the following expression for the total energy, 

£ = |(dr)(r»») = i|(dr)(jV»-V!)G(.,,..': 



i/,dr)/^2.=^(r,r,, 



(16) 



if we integrate by parts and omit the surface term. Integrating over the Green's 
functions in the three regions, given by ( iTal i. ( fTbl i. and ( iTcl i. we obtain for A = A', 
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^^ 1 r. 2 1 ^ Th 3 l+2/(y + A«) 

ASnVJo ^^ l+y/lAa) 96nVjQ ^^ (^/(Afl) + 1)V- 1' ^ ^ 

where the first term is regarded as an irrelevant constant (A is constant so the a can 
be scaled out), and the second term coincides with the massless limit of the energy 
first found by Bordag et al. Il39i , and given in 14411491 . When differentiated with 
respect to a, iT7\ . with A fixed, yields the pressure (fT4b . (We will see below that the 
divergent constant describe the self-energies of the two plates.) 

If, however, we integrate the interior and exterior energy density directly, one 
gets a different result. The origin of this discrepancy with the naive energy is the 
existence of a surface contribution to the energy. To see this, we must include the 
potential in the stress tensor, 

T^'' = ^^',j>^',p - ^-g^' (^d^,pd,,j> + V<p^) , (18) 

and then, using the equation of motion, it is immediate to see that the energy density 
is 

roo ^ ^^o^^o^ _ \Hd'f^ + ^V . (0V0), (19) 

so, because the first two terms here yield the last form in (fTSI l. we conclude that 
there is an additional contribution to the energy. 



i/dS-VG(x,x') 
1 r dco r (dk) ^ d 



I (■ (dk) ^ d 



(20a) 



(20b) 



where the derivative is taken at the boundaries (here x = 0, a) in the sense of the 
outward normal from the region in question. When this surface term is taken into 
account the extra terms incorporated in ( [TtI i are supplied. The integrated formula 
(fT6b automatically builds in this surface contribution, as the implicit surface term in 
the integration by parts. That is, 

E= f{dr){T^°)+E. (21) 

(These terms are slightly unfamiliar because they do not arise in cases of Neumann 
or Dirichlet boundary conditions.) See Fulling ll50l for further discussion. That the 
surface energy of an interface arises from the volume energy of a smoothed interface 
is demonstrated in ||45l, and elaborated in Sect. 12.21 
In the limit of strong coupling, we obtain 

n^ 

limcf^ J, (22) 

A^oo 1440a3' 
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which is exactly one-half the energy found by Casimir for perfectly conducting 
plates m. Evidently, in this case, the TE modes (calculated here) and the TM modes 
(calculated in the following subsection) give equal contributions. 



2.1 TM Modes 



To verify this last claim, we solve a similar problem with boundary conditions that 
the derivative of g is continuous at x = and a. 



<5 / /^ 

ax 



IS contmuous, 



(23a) 



but the function itself is discontinuous. 



^,(x,x) 



ax 



(23b) 



and similarly at jc = 0. (Here the coupling A has dimensions of length.) These bound- 
ary conditions reduce, in the limit of strong coupling, to Neumann boundary condi- 
tions on the planes, appropriate to electromagnetic TM modes: 



A 



dx' 



](x,x') 



= 0. 



(23c) 



=0,a 



It is completely straightforward to work out the reduced Green's function in this 
case. When both points are between the planes, < x,x' <a. 



,(^,y) = i_e-4v-.'l , 1 



Xk 



2cosh)c(x — x') 



Xk / Xk 

IH 

2 I 2 



2kA \ \ 2 

^Kix+J) _j_g-K-(.r+.ir'-2a) 



(24a) 



while if both points are outside the planes, a < x,x', 



Kx,x') = ^e-'^l-^'l 
1 Xk 



-K{x+x'-2a) 



2kA 2 
where the denominator is 



!-■ 



Xk 



Xk^ 
1 + -Ie^ 



(24b) 



Xk 
1 + -I e- 



iKa 



Xk 



(25) 
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It is easy to check that in the strong-coupling limit, the appropriate Neumann 
boundary condition (|23ct is recovered. For example, in the interior region, < 
x^y! < a, 

, ,, coshK'x<coshK'(x> — a) 

limgx,x' = ^ . ■ ^^ -■ (26) 

A^oo fcsinhfca 

Now we can compute the pressure on the plane by computing the xx component 
of the stress tensor, which is given by ( I12al ). so we find 



I _ 1 

■'^ljc=a— 2i 



K- 

A 



2k /Xk^ ^ 



(27a) 



1 



txx\ , = K, (27b) 



and the flux of momentum deposited in the plane x = a is 

and then by integrating over frequency and transverse momentum we obtain the 
pressure: 

In the limit of weak coupling, this behaves as follows: 



dTM 



15 



P'™---— ^A^ (30) 



which is to be compared with ( Il5at . In strong coupling, on the other hand, it has 
precisely the same limit as the TE contribution, (I15bl ). which confirms the expec- 
tation given at the end of the previous subsection. Graphs of the two functions are 
given in Fig.[T] 

For calibration purposes we give the Casimir pressure in practical units between 
ideal perfectly conducting parallel plates at zero temperature: 

n^ ^ 1.30 mPa 

P^ jflC^-—-, -r. (31) 

240^4 (fl/ljum)4 



2.2 Self -energy of Boundary Layer 

Here we show that the divergent self-energy of a single plate, half the divergent term 
in ([TtT i. can be interpreted as the energy associated with the boundary layer We do 
this in a simple context by considering a scalar field interacting with the background 
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0.000 




-0.020 



Fig. 1 TE and TM Casimir pressures between 5-function planes having strength A and separated 
by a distance a. In each case, the pressure is plotted as a function of the dimensionless coupling, 
ka or A/a, respectively, for TE and TM contributions. 



^int==-2^^' 



where the background field a expands the meaning of the 8 function, 



"' 2 ^-^ ^ 2' 

0, Otherwise, 



with the property that h8 = \. The reduced Green's function satisfies 



.9x2 



+ K^ + Xa{x) 



g{x,x') = 5{x-x'). 



This may be easily solved in the region of the slab, —j <x < ^, 



'<-''- 2^1-'^"-"' 4 



XlicoshK'{x + x') 

(K-'-K-)V''^COShK-'(x-y 



(32) 



(33) 



(34) 



(35) 



Here fc' = VK^^+Xli, and 
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(36) 



This result may also easily be derived from the multiple reflection formulas given 
in 1461 . and agrees with that given by Graham and Olum ISTI. 

Let us proceed here with more generality, and consider the stress tensor with an 
arbitrary conformal term 1521 . 

in d + 2 dimensions, d being the number of transverse dimensions, and ^ is an 
arbitrary parameter, sometimes called the conformal parameter. Applying the corre- 
sponding differential operator to the Green's function ( |35l l. introducing polar coor- 
dinates in the (C,^) plane, with ^ = Kcos 0, k = jcsin 9, and 



(sin^e) 



d+l 



(38) 



we get the following form for the energy density within the slab. 



(ro«) = 



2-d-2^-(d+l)/2 



' dKK'' 



r{{d + 2,)/2) Jq k'A 
-{K'-Kf-e-'^'^KA, -5/2 < .=i- < 5/2. 



A/j [(1 -4-^){l+d)K'^- K-2] cosh2K-'x 



(39) 



We can also calculate the energy density on the other side of the boundary, from the 
Green's function forx,x' < —5/2, 



?{x,x') 



2k- 



.^(x+y+5)^/^ 



sinhK''5 



(40) 



and the corresponding energy density is given by 



nOO\ 



d{l~4^{d+l)/d) 



1 



2d+2n{d+i)/2pi^(^^^2)/2)Jo A 



which vanishes if the conformal value of £, is used. An identical contribution comes 
from the region x > 5/2. 

Integrating (r*'") over all space gives the vacuum energy of the slab 



1 



^slab 



2''+27r(''+i)/2r((£/ + 3)/2)Jo 
, , o sinh k' 

+ (1/1)2— — 



dKK" 



K'A 



(K^-KYK^e-^^d 



(42) 



Note that the conformal term does not contribute to the total energy. If we now take 
the limit 5 — )- and h^ °° so that h5 = 1, we immediately obtain the self-energy 
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of a single 5-function plate: 

which for d = 2 precisely coincides with one-half the constant term in (ITtI i. 

There is no surface term in the total Casimir energy as long as the slab is of finite 
width, because we may easily check that -^gl , is continuous at the boundaries 

±-|. However, if we only consider the energy internal to the slab we encounter not 
only the integrated energy density but a surface term from the integration by parts — 
see (l2Tl i. It is the complement of this boundary term that gives rise to Eg, (1431 ). in 
this way of proceeding. That is, as 5 — > 0, 

-/ («fr)/fl'CC'^(r,r)=0, (44) 

J slab J 

SO 

^5=^l.v=-5/2+^l.=5/2' (45) 

with the normal defining the surface energies pointing into the slab. This means that 
in this limit, the slab and surface energies coincide. 

Further insight is provided by examining the local energy density. In this we 
follow the work of Graham and Olum lISTl l53l . From ( [39] l we can calculate the 
behavior of the energy density as the boundary is approached from the inside: 

(7^00. rid+l)Xh 1-4^(^+1)/^ ,,5/2 (46) 

^' ) 2d+^Ki^M)/2r{{d+3)/2) {5-2\x\y ' i-^i^''/^- ^^^^ 

For c/ = 2 for example, this agrees with the result found in fSTJI for ^ = 0: 

Note that, as we expect, this surface divergence vanishes for the conformal stress 
tensor ||52| , where ^ = d/4{d+l). (There will be subleading divergences if d > 2.) 
The divergent term in the local energy density from the outside, dTTl i. as x -^ —5/2, 
is just the negative of that found in ( |46] |. This is why, when the total energy is 
computed by integrating the energy density, it is finite for d <2, and independent 
of ^ . The divergence encountered for d ~2 may be handled by renormalization of 
the interaction potential [l5ll|. 

Note, further, that for a thin slab, close to the exterior but such that the slab still 
appears thin, x^ 5, the sum of the exterior and interior energy density divergences 
combine to give the energy density outside a 5-function potential: 



"•-^spc-"?) 



(*- 5/2)2 (,, + 5/2)2 



' '-^« (48) 
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for small x. Although this limit might be criticized as illegitimate, this result is 
correct for a 5-function potential, and we will see that this divergence structure 
occurs also in spherical and cylindrical geometries, so that it is a universal surface 
divergence without physical significance, barring gravity. 
For further discussion on surface divergences, see Sect. [3] 



3 Surface and Volume Divergences 

It is well known as we have just seen that in general the Casimir energy density 
diverges in the neighborhood of a surface. For flat surfaces and conformal theories 
(such as the conformal scalar theory considered above 1361 . or electromagnetism) 
those divergences are not presentjj In particular. Brown and Maclay 1581 calculated 
the local stress tensor for two ideal plates separated by a distance a along the z axis, 
with the result for a conformal scalar 

{Tn--j^,[^t^t^^gn- (49) 



This result was given more recent rederivations in 15911361 . Dowker and Kennedy 
l60ll and Deutsch and Candelas ||6ri considered the local stress tensor between 
planes inclined at an angle a, with the result, in cylindrical coordinates (f ,r, 0,z), 



(rn = -/^"^ 



llQn^r^ 



0-100 

3 

\0 0-1/ 



(50) 



where for a conformal scalar, with Dirichlet boundary conditions, 

and for electromagnetism, with perfect conductor boundary conditions. 

For a —> we recover the pressures and energies for parallel plates, (I15bl ) and ( |3T1 ). 
(These results were later discussed in ll62l .) 

Although for perfectly conducting flat surfaces, the energy density is finite, for 
electromagnetism the individual electric and magnetic fields have divergent RMS 



' In general, this need not be the case. For example, Romeo and Saharian |54| show that with mixed 
boundary conditions the surface divergences need not vanish for parallel plates. For additional 
work on local effects with mixed (Robin) boundary conditions, applied to spheres and cylinders, 
and corresponding global effects, see I55ll56|[57ll50l . See also Sect. l2.2l and I51II53I . 
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values, 

(E^) ^ -{B^) ^ 1^^ e^O, (53) 

a distance e above a conducting surface. However, if the surface is a dielectric, 
characterized by a plasma dispersion relation, these divergences are softened 

(£2)^J^, -(b2)^1 e^O, (54) 



so that the energy density also diverges 



nOO\ 



1 



(r-)^_, £^0. (55) 

The null energy condition {n^^n^ = 0) 

T^^rif^riy > (56) 

is satisfied, so that gravity still focuses light. 

Graham 1651 166| examined the general relativistic energy conditions required 
by causality. In the neighborhood of a smooth domain wall, given by a hyperbolic 
tangent, the energy density is always negative at large enough distances. Thus the 
weak energy condition is violated, as is the null energy condition ( |56] ). However, 
when ( l56b is integrated over a complete geodesic, positivity is satisfied. It is not 
clear if this last condition, the Averaged Null Energy Condition, is always obeyed in 
flat space. Certainly it is violated in curved space, but the effects always seem small, 
so that exotic effects such as time travel are prohibited. 

However, as Deutsch and Candelas 11611 showed many years ago, in the neigh- 
borhood of a curved surface for conformally invariant theories, (T^v) diverges as 
e^^, where e is the distance from the surface, with a coefficient proportional to the 
sum of the principal curvatures of the surface. In particular they obtain the result, in 
the vicinity of the surface, 

(7),v) ^ e-'4l^ + e-^4^J + E-'4l\ (57) 

and obtain explicit expressions for the coefficient tensors T^y and Tr^y in terms of 
the extrinsic curvature of the boundary. 

For example, for the case of a sphere, the leading surface divergence has the 
form, for conformal fields, for r = fl + e, e— >0 

\ 

(58) 

in spherical polar coordinates, where the constant is A = l/7207r^ for a scalar field 
satisfying Dirichlet boundary conditions, or A = 1 /607r^ for the electromagnetic 







/2/fl 


7mv> = 


A 



fl 
I flsin^ 
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field satisfying perfect conductor boundary conditions. Note that (ISST i is properly 
traceless. The cubic divergence in the energy density near the surface translates 
into the quadratic divergence in the energy found for a conducting ball 1671. The 
corresponding quadratic divergence in the stress corresponds to the absence of the 
cubic divergence in (7)-^). 

This is all completely sensible. However, in their paper Deutsch and Candelas 
ll6n expressed a certain skepticism about the validity of the result of 1681 for the 
spherical shell case (described in part in Sect. 14.2b where the divergences cancel. 
That skepticism was reinforced in a later paper by Candelas ||69l, who criticized the 
authors of [681 for omitting 5 function terms, and constants in the energy. These ob- 
jections seem utterly without merit. In a later critical paper by the same author ITOl . 
it was asserted that errors were made, rather than a conscious removal of unphysical 
divergences. 

Of course, surface curvature divergences are present. As Candelas noted 1691 iTOl . 
they have the form 

E=E^ fdS + E^ [dS{Ki + K-2) +Ef [dS{Ki - K2f+Efj f dSKiK2 + ... (59) 

where iCi and K2 are the principal curvatures of the surface. The question is to what 
extent are they observable. After all, as has been shown in 13811361 and in Sect. 12.21 
we can drastically change the local structure of the vacuum expectation value of the 
energy-momentum tensor in the neighborhood of flat plates by merely exploiting the 
ambiguity in the definition of that tensor, yet each yields the same finite, observable 
(and observed!) energy of interaction between the plates. For curved boundaries, 
much the same is true. A priori, we do not know which energy-momentum ten- 
sor to employ, and the local vacuum-fluctuation energy density is to a large extent 
meaningless. It is the global energy, or the force between distinct bodies, that has 
an unambiguous value. It is the belief of the author that divergences in the energy 
which go Uke a power of the cutoff are probably unobservable, being subsumed in 
the properties of matter Moreover, the coefficients of the divergent terms depend on 
the regularization scheme. Logarithmic divergences, of course, are of another class 
l40l . Dramatic cancellations of these curvature terms can occur. It might be thought 
that the reason a finite result was found for the Casimir energy of a perfectly con- 
ducting spherical shell 1131120 68 1 is that the term involving the squared difference 
of curvatures in ( i59b is zero only in that case. However, it has been shown that at 
least for the case of electromagnetism the corresponding term is not present (or has 
a vanishing coefficient) for an arbitrary smooth cavity 1211 . and so the Casimir en- 
ergy for a perfectly conducting ellipsoid of revolution, for example, is finite^ This 
finiteness of the Casimir energy (usually referred to as the vanishing of the sec- 
ond heat-kernel coefficient fTTl ) for an ideal smooth closed surface was anticipated 
already in l20l . but contradicted by 1611 . More specifically, although odd curva- 



^ The first steps have been made for calculating the Casimir energy for an ellipsoidal boundary 
1341 1351 . but only for scalar fields since the vector Helmholtz equation is not separable in the 
exterior region. 
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ture terms cancel inside and outside for any thin shell, it would be anticipated that 
the squared-curvature term, which is present as a surface divergence in the energy 
density, would be reflected as an unremovable divergence in the energy. For a closed 
surface the last term in (|59] l is a topological invariant, so gives an irrelevant constant, 
while no term of the type of the penultimate term can appear due to the structure of 
the traced cylinder expansion ll50l . 



4 Casimir Forces on Spheres via 5-Function Potentials 

This section is an adaptation and an extension of calculations presented in 1451 |46l . 
This investigation was carried out in response to the program of the MIT group 
Il4ni42ll43ll44l|49l . They first rediscovered irremovable divergences in the Casimir 
energy for a circle in 2 + 1 dimensions first discovered by Sen 1721 [73ll . but then 
found divergences in the case of a spherical surface, thereby casting doubt on the 
vahdity of the Boyer calculation ITSl . Some of their results, as we shall see, are 
spurious, and the rest are well known l 40| . However, their work has been valuable 
in sparking new investigations of the problems of surface energies and divergences. 
We now carry out the calculation we presented in Sect. |2|in three spatial dimen- 
sions, with a radially symmetric background 

^M^-~8{r-a)^\x), (60) 

which would correspond to a Dirichlet shell in the limit X ^> °°. The scaling of the 
coupling, which here has dimensions of length, is demanded by the requirement 
that the spatial integral of the potential be independent of a. The time-Fourier trans- 
formed Green's function satisfies the equation [k^ ~ — o^) 

-V^ + K^ + ^5{r-a) ^(r,r') = 5(r-r'). (61) 

fl^ J 

We write ^ in terms of a reduced Green's function 

^(r,r') = '£g,{r,r')YUQ)Y,:„iQ^), (62) 

Im 



where gi satisfies 



1 d 2d , /(/ + 1) , .^2 , ^ ,r , 

- — — r — + , — + K +—5{r-a) 

n dr dr r" a^ 



g,{r,r') = ^5{r~r'). (63) 



We solve this in terms of modified Bessel functions, /y (x), Ky (x), where v = / + 1 /2, 
which satisfy the Wronskian condition 

/; {x)Ky (x) - K'y {x)Iy {x) = -. (64) 

X 
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The solution to ( |63T l is obtained by requiring continuity of gi at each singularity, at 
r' and a, and the appropriate discontinuity of the derivative. Inside the sphere we 
then find (0 < r,r' < a) 



gl{r,r') 



Krf 



ei{Kr>)si{Kr<) :^si{Kr)s,{Kr') 



ej{Ka) 



^^" l + -^.9,((Cfl)e,((Cfl) 



(65) 



Here we have introduced the modified Riccati-Bessel functions. 



■s/W = \Hrh+\l2{x), ei(x) = \ —K,+xi2{x) 



(66) 



Note that d&Sl l reduces to the expected Dirichlet result, vanishing as r — > a, in the 
limit of strong coupling: 



\\m gi{r,r') = —-I 
A^oo Krr 



ei{Ka) , . , i' 
e,{Kr>)s,{Kr<) '- — -si{Kr)si{Kr') 



si{Ka) 
When both points are outside the sphere, r, r' > a, we obtain a similar result: 



8l{r,r')- 



Krr' 



A 
e,(K-r>>,(K-r<) - -—^ei{Kr)e,{Kr') 



sj{Ka) 



l + -^si{Ka)ei{Ka) 



(67) 



(68) 



which similarly reduces to the expected result as A ^ oo. 

Now we want to get the radial-radial component of the stress tensor to extract the 
pressure on the sphere, which is obtained by applying the operator 



d,d,-\i-d'd'' + V.V')^i 



drd,j -K" - 



Ki + i) 



(69) 



to the Green's function, where in the last term we have averaged over the surface of 
the sphere. Alternatively, we could notice that ITMI 



V-V'P/(cos7) 



1(1 + 1) 



(70) 



y-^O 



where 7 is the angle between the two directions. In this way we find, from the 
discontinuity of {T,-,-) across the r = a surface, the net stress 



.y 



-3L(2' + i)/ ^ — 



2na 



1=0 







1 



Xaei{x)si{x) 



(71) 



(Notice that there was an error in the sign of the stress, and of the scaling of the 
coupling, in B3]|46l .) 

The same result can be deduced by computing the total energy ( fTSI l. The free 
Green's function, the first term in (l65T l or (l68T l. evidently makes no significant con- 
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tribution to the energy, for it gives a term independent of the radius of the sphere, a, 
so we omit it. The remaining radial integrals are simply 

fdysjiy) = ^ [{x^ + lil + l))sj{x)+xs,ix)s',ix)~x^s'^{x)] , (72a) 

J^yeKy) = -^ [{x^ + l{l + l))ejix)+xei{x)e',{x)-xy^ix)] . (72b) 



(73) 



Then using the Wronskian (|64] |. we find that the Casimir energy is 



E = y(2l + l)[ djcx — In 



1 + -Iv{x)Ky{x) 



If we differentiate with respect to a we immediately recover the force dTTl i. This 
expression, upon integration by parts, coincides with that given by Barton llTSl . and 
was first analyzed in detail by Scandurra 1761 . This result has also been rederived 
using the multiple-scattering formalism 1251 . For strong coupling, it reduces to the 
well-known expression for the Casimir energy of a massless scalar field inside and 
outside a sphere upon which Dirichlet boundary conditions are imposed, that is, that 
the field must vanish at r = a: 

Hm £ = - J- y (2/ + 1 ) / djcx -^ In \Iy (x)Ky (x)] , (74) 

l^oo 2na f^^ Jo ax 

because multiplying the argument of the logarithm by a power of x is without effect, 
corresponding to a contact term. Details of the evaluation of (|74] | are given in l36l . 
and will be considered in Sect. g^below. (See also I77ll78ll79l .) 

The opposite limit is of interest here. The expansion of the logarithm is immedi- 
ate for small /L. The first term, of order A, is evidently divergent, but irrelevant, since 
that may be removed by renormalization of the tadpole graph. In contradistinction 
to the claim of ||42l|43]|44l|49l, the order X^ term is finite, as estabHshed in [361. 
That term is 

£^^'^ = ^E(2/ + l)j[^"ckx^[//+i/2W^;+i/2W]'- (75) 



The sum on / can be carried out using a trick due to Klich l80l : The sum rule 

oo 

Y,i2l + l)e,ix)s,{y)P,icose) = ^e^P, (76) 

1=0 P 

where p = ^x^ +y^ — 2xycos0, is squared, and then integrated over 9, according 

to 

"1 2 

d(cos0)P,(cos0)P//(cos0) = 5/// . (77) 

1 2/ + 1 

In this way we learn that 
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2 /■4^dw 



l^i2l + l)ej{x)s^ix)^- -e-- 



1=0 



2 Jo w 



(78) 



Although this integral is divergent, because we did not integrate by parts in dTSl l. 
that divergence does not contribute: 



.(A^ 



d.i.Ar^e-"'= 



4;rfl- Jo 2 djc Jo w 



32;rfl3 ' 



(79) 



which is exactly the result (4.25) of | 

However, before we wax too euphoric, we recognize that the order A^ term ap- 
pears logarithmically divergent, just as ll44ll and 1491 claim. This does not signal a 
breakdown in perturbation theory. Suppose we subtract off the two leading terms. 



— f;(2/+i) Hducx— 

Ina f^Q Jo dx 

^ 32na^' 



In 1 



-IvK\, 



X X 

-aIyKy + -^{lyKyy 

a 2a^ 



(80) 



To study the behavior of the sum for large values of /, we can use the uniform 
asymptotic expansion (Debye expansion), for v — > 0°, 



/v(x) 



27rv 



^Ml+I 



Uk{t) 



^vw-/|e--(i+i:(-i)^^~ 



(81) 



where 



x = mz, f = l/VTT?, Tj(z) = \/l+z- + ln 



l+W- 



d77 _ 1 

dz zt' 



The polynomials in t appearing in dSTT i are generated by 



(82) 



2 8 Jo 

We now insert these expansions into (ISOl l and expand not in A but in v; the leading 
term is 

,3, A^ ^ 1 r dz x^ 



£(A^') 



•,4L ,, 



247rfl\t'o ^ -^0 (l+z2)3/2 24na 



4^(1) 



(84) 



Although the frequency integral is finite, the angular momentum sum is divergent. 
The appearance here of the divergent i^(l) seems to signal an insuperable barrier 
to extraction of a finite Casimir energy for finite A . The situation is different in the 
Umit A -> 00 — See Sect. 1421 
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This divergence has been known for many years, and was first calculated explic- 
itly in 1998 by Bordag et al. ||40| . where the second heat kernel coefficient gave an 
equivalent result, 

4-671 a^ s 
A possible way of dealing with this divergence was advocated in (|26|. More re- 
cently, Bordag and Vassilevich lISTI have reanalyzed such problems from the heat 
kernel approach. They show that this 0{X^) divergence corresponds to a surface ten- 
sion counterterm, an idea proposed by me in 1980 ll82l[83]| in connection with the 
zero-point energy contribution to the bag model. Such a surface term corresponds 
to X fixed, which then necessarily implies a divergence of order A^. Bordag argues 
that it is perfectly appropriate to insert a surface tension counterterm so that this 
divergence may be rendered finite by renormalization. 



4.1 TM Spherical Potential 

Of course, the scalar model considered in the previous subsection is merely a toy 
model, and something analogous to electrodynamics is of far more physical rele- 
vance. There are good reasons for believing that cancellations occur in general be- 
tween TE (Dirichlet) and TM (Robin) modes. Certainly they do occur in the classic 
Boyer energy of a perfectly conducting spherical shell ifTSl |20l |68l . and the indi- 
cations are that such cancellations occur even with imperfect boundary conditions 
ITTSI . Following the latter reference, let us consider the potential 

=S^nt-iA-|-5(r-fl)(/)2(x). (86) 

2 r or 

Here A again has dimensions of length. In the limit A ^- oo this corresponds to TM 
boundary conditions. The reduced Green's function is thus taken to satisfy 



r^ or or r^ r or 



gl{r,r')^\5{r-r'). (87) 



At r = r' we have the usual boundary conditions, that gi be continuous, but that its 
derivative be discontinuous. 



id 

r -^gi 
ar 



r=r'+ 

= -1, (88) 

r=r'- 



while at the surface of the sphere the derivative is continuous. 



d 



■■ 0, (89a) 
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r=a-\- 



X d 



-rgi 



(89b) 



Equations ( |89at and ( |89b| ) are the analogues of the boundary conditions ( I23ab . ( I23bl i 
treated in Sect. 12.11 

It is then easy to find the Green's function. When both points are inside the 
sphere, 



r,r <a: gi{r,r) ■ 



1 
Krr' 



■ XK[e',{Ka)Ysi{Kr)si{Kr') 

si{Kr^)ei{Kry) ; — — t- — r — 

^ ^^ ' l+XKe'i{Ka)s'i{Ka) 



and when both points are outside the sphere. 



r,r'>a: gi{r,r') 



Krr 



s,{Kr<)ei{Kr>) ^^ 



XK[s'i{Ka)Yei{Kr)ei{Kr') 
Ke\{Ka)s'i{Ka) 



(90a) 



(90b) 



It is immediate that these supply the appropriate Robin boundary conditions in the 
X ^ °° limit: 

\im—rgi =0. (91) 

X^odr ,.^„ 

The Casimir energy may be readily obtained from ( fT6] l, and we find, using the 
integrals ( I72al ). ( I72bb 



1 °° r°° d 

Vill + l) / djcx— In 

'~' Jo djc 



iKa 



1=0 



X 
1 H — xe'i{x)s'i{x) 



(92) 



The stress may be obtained from this by applying —d/da, and regarding X as con- 
stant, or directly, from the Green's function by applying the operator. 



_ 1 

2l 






gl 



(93) 



which is the same as that in ( |69l ). except that 



r 



(94) 



appropriate to TM boundary conditions (see 11841 . for example). Either way, the total 
stress on the sphere is 



y : 



' ti2i+i)rdxx- [^'W^'Wi' 



iKO^ 



>^..Jl^\.ll 



/=o -/o \ + ^xe',{x)s\{x) 



(95) 



The result for the energy ( |92] i is similar, but not identical, to that given by Barton 

ca. 
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Suppose we now combine the TE and TM Casimir eneraies, (TTSll and (|92li: 


£™+£TM_ 1 f (2/ + 1) Tdjcx'^ln 
2na /■^Q^ Jo dx 


\ a X J \ a _ 




(96) 


In the limit X ^- oo this reduces to the familiar expression for the perfectly conduct- 


ing spherical shell |68|: 


lim£ = y(2l + l)[ dxx 


\ei e, s, sj 


(97) 



Here we have, as appropriate to the electrodynamic situation, omitted the / = 
mode. This expression yields a finite Casimir energy, as we will see in Sect. 14.21 
What about finite A? In general, it appears that there is no chance that the divergence 
found in the previous section in order A ^ can be cancelled. But suppose the coupling 
for the TE and TM modes are different. If A^^A^'^ = 4fl^, a cancellation appears 
possible, as discussed in 



4.2 Evaluation of Casimir Energy for a Dirichlet Spherical Shell 

In this subsection we will evaluate the above expression (|74] | for the Casimir energy 
for a massless scalar in three space dimensions, with a spherical boundary on which 
the field vanishes. This corresponds to the TE modes for the electrodynamic situ- 
ation first solved by Boyer ifTsl l20l l68l . The purpose of this section (adapted from 
ll36ll46l ) is to emphasize anew that, contrary to the implication of Ii42n43, , 44,. 49l . 
the corresponding Casimir energy is also finite for this configuration. 

The general calculation in D spatial dimensions was given in llTTl : the pressure 
is given by the formula 



^ 2/+£)-2 r /+D-2 r d r , , , , 7 m 
P= y -^^ , ,; ^ „ , '- / dxx— In \ly (x)Kv (x)x^-'^] 

Here v = / — 1 + D/2. For D = 3 this expression reduces to 



(98) 



P = -^^ £(2/ + l)^"d«Ain [/,+i/2W/:,+i/2(x)/x] . (99) 

This precisely corresponds to the strong limit X ^ °° given in (|74l l. if we recall the 
comment made about contact terms there. In iTTl we evaluated expression (|98] | by 
continuing in D from a region where both the sum and integrals existed. In that way, 
a completely finite result was found for all positive D not equal to an even integer. 

Here we will adopt a perhaps more physical approach, that of allowing the time- 
coordinates in the underlying Green's function to approach each other, temporal 
point-sphtting, as described in ll68l . That is, we recognize that the x integration 
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above is actually a (dimensionless) imaginary frequency integral, and therefore we 
should replace 

/■OO 1 /"OO 

J^ dxf{x) = -J Jye'y^f{\y\), (100) 

where at the end we are to take 5^-0. Immediately, we can replace the x^^ inside 
the logarithm in ( |99l ) by x, which makes the integrals converge, because the differ- 
ence is proportional to a 5 function in the time separation, a contact term without 
physical significance. 

To proceed, we use the uniform asymptotic expansions for the modified Bessel 
functions, dsTl ). This is an expansion in inverse powers of v = / + 1/2, low terms in 
which turn out to be remarkably accurate even for modest /. The leading terms in 
this expansion are, using dSTT i. 

ln[x/,+i/2(x)/:/+i/2(x)] ~ln| + ^§(0 + :^MO + ---, (101) 



g(t) = ^{t^-6t^ + 5t^), (102a) 

h{t) = ;^(13f'*-284f''+1062f**-1356f"' + 565f^^). (102b) 

64 

The leading term in the pressure is therefore 

(103) 
where in the last step we have used the formal zeta function evaluatioqj 

OO 

^v-^ = (2^-l)CW. (104) 

1=0 

Here the rigorous way to argue is to recall the presence of the point-splitting factor 
givzS ^jj(j jQ carry out the sum on / using 



1=0 

so 



:e--^^ = -i-^, (105) 

2i smzo/2 



fvV^-^5 = ^ I = L( ? + —^ ] (106) 

,t-o d(z5)2 2sinz5/2 8^ smh5/2 sinzd/lj ' ^ ' 



' Note that the corresponding TE contribution the electromagnetic Casimir pressure would not be 
zero, for there the sum starts from / = 1 . 
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Then Pq is given by the divergent expression 

which we argue is zero because the integrand is odd, as justified by averaging over 
contours passing above and below the pole at z = 0. 

The next term in the uniform asymptotic expansion (llQlt . that involving g, like- 
wise gives zero pressure, as intimated by ( |104t . which vanishes at s = 0. The same 
conclusion follows from point splitting, using ( llQSb and arguing that the resulting 
integrand ~' z^t^g'{t)/z8 is odd in z- Again, this cancellation does not occur in the 
electromagnetic case because there the sum starts at Z = 1 . 

So here the leading term which survives is that of order V^^ in dlOlt . namely 

^2 = T^ E \ rdz/i(0, (108) 

where we have now dropped the point-splitting factor because this expression is 
completely convergent. The integral over z is 

dzh(t) = — (109) 

^ ' 32768 

and the sum over Hs 3i^(2) = 7:^/2, so the leading contribution to the stress on the 
sphere is 

.. . 2 357r2 0.00527094 

65536fl- a^ 

Numerically this is a terrible approximation. 

What we must do now is return to the full expression and add and subtract the 
leading asymptotic terms. This gives 

^ = J^2-^£(2/ + l)i?/, (Ill) 



where 

R, = Qi+ I dx • 

Jo 



Inzt + -jg{t) + -jh{t) 



(112) 



where the integral 

Q, = - [ dx\n[2xh{x)Kv{x)] (113) 

Jo 

was given the asymptotic form in ||77][38l (/ >• 1): 

VTT n 35n 565;r 1208767?: 



Qi 



2 128v 32768v3 1048577v5 2147483648v^ 
138008357;: 



137438953472v'^ 



(114) 
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The first two terms in ( I114l i cancel the second and third terms in ( 1112b . of course. 
The third term in ( I114l i corresponds to h{t), so the last three terms displayed in ( I114l i 
give the asymptotic behavior of the remainder, which we call w(v). Then we have, 
approximately. 



^«^2--^i:v/;,--^ £ VW(V). 



%a- f-'^ Tta- ,^„^j 



(115) 



For n — I this gives S^ « 0.002 852 78/fl^, and for larger n this rapidly approaches 
the value first given in fTTl . and rederived in 1178117911851 

^™ = 0.0028 17/fl2, (116) 

a value much smaller than the famous electromagnetic result | 

^- = ^:^, (117) 

because of the cancellation of the leading terms noted above. Indeed, the TM con- 
tribution was calculated separately in 1841 . with the result 

J5^™ = -0.02204-r, (118) 

and then subtracting the I = modes from both contributions we obtain ( 11171 ). 

48fl^ a'^ 



4.3 Surface Divergences in the Energy Density 

The following discussion is based on 1741. Using ( fTOl i. we immediately find the 
following expression for the energy density inside or outside the sphere; 



,oon / 'dic^2/+l 



^")- .i^L 



2a: /■^^ Ak 



.2, ,,. , /(^ + 1) 



-K^ + d,-d,j 



8l(.r,r') 



-2^\^r^^g,{rA, (120) 

n or or J 

where ^ is the conformal parameter as seen in ( l37l i. To find the energy density in 
either region we insert the appropriate Green's functions (i65i l or ( i68l l. but delete the 
free part, 

g'l = -^^ii^'-<>i{^'->). (121) 
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which corresponds to the bulk energy which would be present if either medium filled 
all of space, leaving us with for r > a 



«(r) = -(l-4^; 



■dK-^ 2Z + 1 






'o 2nf-'^ An \ + ^e,{Ka)si{Ka)\. Kr 



ej{Kr) 



-K' 



1(1+1) 



're,{Kr)e',{Kr) + ^e'^{Kr) 



1+4^ 
1-4^ 

(122) 



Inside the shell, r < a, the energy is given by a similar expression obtained from 
(1122b by interchanging e/ and si . 

We want to examine the singularity structure as r — > a from the outside. For this 
purpose we use the leading uniform asymptotic expansion, I -^ °o, obtained from 
(EB 



ei{x) - y^e "'', si{x) - 2^^"''' 



e'lix) 



-vr] 



Zt 



s',{x)' 



1 1 

2^' 



,vn 



(123) 



where V = / + 1/2, andz, t, and Tj are given in ( |82] |. The coefficient of e/(KT)e/(KT') 
occurring in the 5-function potential Green's function ( |68] ), in strong and weak cou- 
pUng, becomes 



A 
a 



0: 



si{Ka) 
ei{Ka)' 



(124a) 
(124b) 



In either case, we carry out the asymptotic sum over angular momentum using 
(1123b and the analytic continuation of (11051 ) 



Here (r sa a) 



X 





1=0 


1 






2sinhf ' 




'niz) 


-"im 


dz r 


2r~a 
t r 



(125) 



(126) 



The remaining integrals over z are elementary, and in this way we find that the 
leading divergences in the energy density are as r ^ «+, 



X 
a 



A") 



1 l-6.g 

167r2(r-fl)4' 

A\" r(4-n) 

a 967r2fl4 ^ ^' 



(127a) 



4-n 



, «<4,(127b) 
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where the latter is the leading divergence in order n. These results clearly seem to 
demonstrate the virtue of the conformal value of ^ = 1 /6; but see below. (The value 
for the Dirichlet sphere ( I127al ) first appeared in |i6l|; it more recently was rederived 
in lISTl . where, however, the subdominant term, the leading term if (^ = 1/6, namely 
(I130l l, was not calculated. Of course, this result is the same as the surface divergence 
encountered for parallel Dirichlet plates ll36l[38l .) The perturbative divergence for 
w = 1 in (I127bb is exactly that found for a plate — see ( |48] |. 

Thus, for ^ = 1/6 we must keep subleading terms. This includes keeping the 
subdominant term in ;if O 

2 r — a f r — a\ 
,.-^ + r(--j , (128) 

the distinction between t{z) and I = t{z = za/r), 

t r — a 

zt^zt-th , (129) 

r 

as well as the next term in the uniform asymptotic expansion of the Bessel functions 
(ISTI ). Including all this, it is straightforward to recover the well-known result ( fSSl l 
ll6n for strong coupling (Dirichlet boundary conditions): 

a 360;r^ fl(r— fl)^ 

Following the same process for weak coupling, we find that the leading divergence 
in order «, 1 < « < 3, is (r — !• azt) 

^-0: «<"^-(^)"T^^(^(«-l)(«+2)r(3-„). (131) 

Note that the subleading <9(A) term again vanishes. Both Eqs. ( 11301 ) and ( 1131b apply 
for the conformal value ^ = 1/6. 



4.4 Total Energy and Renormalization 

As discussed in 1741 we may consider the potential, in the spirit of 

^M^~^<P^o{r), (132a) 

where 



Note there is a sign en'or in (4.8) of 1741 . 
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{0, r<a^, 
h,a^ <r<a+, (132b) 

0, fl+ < r. 

Here a± = a± 5/2, and we set h5 = 1. That is, we have expanded the 5-function 
shell so that it has finite thickness. 

In particular, the integrated local energy density inside, outside, and within the 
shell is Ein, Eom, and E^h, respectively. The total energy of a given region is the sum 
of the integrated local energy and the surface energy ( I20ab bounding that region 
(^ = 1/6): 

Ein=Ei„+E^, (133a) 

^out = £'out +E+, (133b) 

£sh=£sh +£;+£-, (133c) 

where E± is the outside (inside) surface energy on the surface at r = a±, while E'^ is 
the inside (outside) surface energy on the same surfaces. E-^, Eom, and E^h represent 
/(dr)(r''") in each region. Because for a nonsingular potential the surface energies 
cancel across each boundary. 



-£+=0, E^+E'_=0, (134) 



the total energy is 



■^sh — ^in I ^out I ^sh- 



(135) 



In the singular thin shell limit, the integrated local shell energy is the total surface 
energy of a thin Dirichlet shell: 

£sh =£++£- 7^0. (136) 



See the remark at the end of Sect. 12. 21 This shell energy, for the conformally coupled 
theory, is finite in second order in the coupling (in at least two plausible regulariza- 
tion schemes), but diverges in third order We showed in ll74l that the latter precisely 
corresponds to the known divergence of the total energy in this order. Thus we have 
established the suspected correspondence between surface divergences and diver- 
gences in the total energy, which has nothing to do with divergences in the local 
energy density as the surface is approached. This precise correspondence should 
enable us to absorb such global divergences in a renormalization of the surface en- 
ergy, and should lead to further advances of our understanding of quantum vacuum 
effects. We will elaborate on this point in the following. 
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5 Semitransparent Cylinder 

This section is based on |!37l. We consider a massless scalar field in a 5-cylinder 
background, 

Xnt^ - ^5 ir-a)(^\ (137) 

a being the radius of the "semitransparent" cylinder The massive case was earlier 
considered by Scandurra ISSll . We will continue to assume that the dimensionless 
coupling A > to avoid the appearance of negative eigenfrequencies. The time- 
Fourier transform of the Green's function satisfies 

[-V2-a)2 + A5(r-fl)]^(r,r') = 5(r-r'). (138) 

Adopting cylindrical coordinates, we write 

(139) 





J In ,„f-„ 


^e""^'^-«>'^gUr,r';k), 


where the reduced Green's function satisfies 




— r— + K-^ + ^ + -5(r-fl) 

r dr ar r^ a 


g„,{r,r';k) = U{r~/) 



(140) 
r- a J r 

where k^ — k^ — aP". Let us immediately make a EucHdean rotation, 

co^iC, (141) 

where i^ is real, so K" is likewise always real. Apart from the 5 functions, this is the 
modified Bessel equation. 

Because of the Wronskian ( l64l i satisfied by the modified Bessel functions, we 
have the general solution to ( 1140b as long as r 7^ a to be 

^,„(r,r';fc) = /,„(fCr<)/i:„,()Cr>) +A(r')/„,(fCr) +Z?(r')/i:„,()Cr), (142) 

where A and B are arbitrary functions of r'. Now we incorporate the effect of the 5 
function at r = a in (1140b . It implies that g,„ must be continuous at r = a, while it 
has a discontinuous derivative. 



—g,n{r,r';k) 
dr 



r=a+ -| 

(fl,r';fc), (143) 



a 



from which we rather immediately deduce the form of the Green's function inside 
and outside the cylinder: 

r,r' <a: g,„ {r,r';k) ^ /,„ ( Kr< )K,n ( Kr^ ) 
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I + Xl,„{Ka)K,„{Ka) 
r,r' >a: g,„ {r,r';k) = I,„ ( KT< )K„, ( Kr> ) 

l+Xl,„{Ka)K,„{Ka) 



I,„iKr)I,„{K/), (144a) 
K,„{Kr)K,„iKr'). (144b) 



Notice that in the limit A — > 0° we recover the Dirichlet cylinder result, that is, that 
gm vanishes at r = a. 



5.1 Cylinder Pressure and Energy 

The easiest way to calculate the total energy is to compute the pressure on the cylin- 
drical walls due to the quantum fluctuations in the field. This may be computed, at 
the one-loop level, from the vacuum expectation value of the stress tensor, 

(145) 
which we have written in a Cartesian coordinate system. Here we have again in- 
cluded the conformal parameter ^, which is equal to 1/6 for the stress tensor that 
makes conformal invariance manifest. The conformal term does not contribute to the 
radial-radial component of the stress tensor, however, because then only transverse 
and time derivatives act on G{x,x), which depends only on r. The discontinuity of 
the expectation value of the radial-radial component of the stress tensor is the pres- 
sure of the cylindrical wall: 

-* ^ \^rr)m \-^rr)out 

1 ^ r^r-^^_ Xk- 



X [Kl{Ka)l';^,{Ka)-ll{Ka)K',l{Ka)\ 

Y, I dk I di;-- — ln[l+Xl,„{Ka)K,„iKa)], (146) 



where we have again used the Wronskian ( |64l i . Regarding ka and l^a as the 
two Cartesian components of a two-dimensional vector, with magnitude x= Ka = 
y^V + ^^fl2, we get the stress on the cylinder per unit length to be 

y = 27taP= t/ dxx^ y — ln\l + XI Jx)Kjx)] , (147) 

4-71 a^ Jo .„_ „djt: 



" m=— o 



which possesses the expected Dirichlet limit as A — )- 0°. The corresponding expres- 
sion for the total Casimir energy per unit length follows by integrating 
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y = -^S, (148) 

da 

that is, 

S = --^rAxx^ £ ^\n[l+Xl,„{x)K,^{x)]. (149) 

This expression, the analog of (|73] | for the spherical case, is, of course, completely 
formal, and will be regulated in various ways, for example, with an analytic or ex- 
ponential regulator as we will see in the following, or by using zeta-function regu- 
larization | f37l . 

Alternatively, we may compute the energy directly from the general formula (fTSI l. 
To evaluate ( fTSI l in this case, we use the standard indefinite integrals over squared 
Bessel functions. When we insert the above construction of the Green's function 
(I144al i and (I144bl i. and perform the integrals over the regions interior and exterior 
to the cylinder we obtain ( |149l l immediately. 



5.2 Weak-coupling Evaluation 

Suppose we regard A as a small parameter, so let us expand (|149t in powers of X. 
The first term is 

'^^'^ = -^ E r^^''^K„,{x)Ux). (150) 

8;rfl^ ,„f^^Jo dx 

The addition theorem for the modified Bessel functions is 

Ko{kP)= £ e™(^-^V„,(fcp)/,„(fcp'), p>p', (151) 

where P = ^/p^ + p'^ — 2pp'cos(0 — ^'). If this is extrapolated to the limit p' = p 
we conclude that the sum of the Bessel functions appearing in ( 11501 ) is Ko{Q), that 
is, a constant, so there is no first-order contribution to the energy. For a rigorous 
derivation of this result, see |[37l . 

We can proceed the same way to evaluate the second-order contribution, 

^^'^ = T^ r^'' J- t lli-)Kl{-)- (152) 

lona^ Jo djc,„r^„ 

By squaring the sum rule ( 11511 ), and taking the limit p' -^ p, we evaluate the sum 
over Bessel functions appearing here as 

£ ll{x)Kl{x)^ r^Ki{2xsm(p/2). (153) 

m=^oo Jo 2% 
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Then changing the order of integration we find that the second-order energy can be 
written as 

^'''=-^2 r-^ fdzz^W, (154) 

where the Bessel-function integral has the value 1/2. However, the integral over (p 
is divergent. We interpret this integral by adopting an analytic regularization based 
on the integral ||3T1 

^2 



^V 2J r 1 + 4) 



which holds for Res > —1. Taking the right-side of this equation to define the (p 
integral for all s, we conclude that the (p integral in (11541) . and hence the second- 
order energy t?'^', is zero. 



5.2.1 Numerical Evaluation 

Given that the above argument evidently formally omits divergent terms, it may be 
more satisfactory, as in [31], to offer a numerical evaluation of (f'^'. (The corre- 
sponding argument for ^''' is given in fJTl.) We can very efficiently do so using the 
uniform asymptotic expansions dSTT l. Thus the asymptotic behavior of the product 
of Bessel functions appearing in (11521 ) is 



^;^W^^W-^(l + L^)- (156) 



The first three polynomials occurring here are 



(2 

ri(0 = -(l-6f2 + 5f4), (157a) 

A 
r2{t) = -— (7-148f2 + 554f'*-708f^ + 295f**), (157b) 

16 

'■3(0 = — (36-1666f2+13775f'^-44272f^ 



16 



-67162r** -48510?'°+ 13475f'^). (157c) 



We now write the second-order energy (I152l i as 
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rk{t) 



m 



2k 



-2 ^ / dxx 



m=lk=l 

2 /„\ 1-2 



/2(x)/:2(x)--^ 1 + ^ 



'■i(0 



4m2 



/t=l 






(158) 



In the final integral z ~ x/m. The successive terms are evaluated as 

A2 



^(2) 



l(r+ln4)-iln2.-M + ZC(l)_iim 
4^^^ '^ 4 48 1920 16128 



87ra2 

+0.000864 + 0.000006 



8;rfl- 



■(0.000000), 



(159) 



where in the last term in ( 11581 ) only the m = 1 and 2 terms are significant. Therefore, 
we have demonstrated numerically that the energy in order A ^ is zero to an accuracy 
of better than 10"''. 

The astute reader will note that we used a standard, but possibly questionable, 
analytic regularization in defining the second term in ([158), where the initial sum 
and integral are only defined for 1 < s < 2, and then the result is continued to s = 0. 
Alternatively, we could follow ifSTI and insert there an exponential regulator in each 
integral of e"^ , with 5 to be taken to zero at the end of the calculation. For m =^ 
X becomes mz, and then the sum on m becomes 






nzS 



1 



e"^ - 1 ' 



Then when we carry out the integral over z we obtain for that term 



n 
85 



1 



■InlTt. 



(160) 



(161) 



Thus we obtain the same finite part as above, but in addition an explicitly divergent 
term 

1 2 

(162) 



^(2) 
°div 



64fl25 ■ 

If we think of the cutoff in terms of a vanishing proper time T, 5 = r/a, this diver- 
gent term is proportional to I /a, so the divergence in the energy goes like L/a, if L 
is the (very large) length of the cylinder This is of the form of the shape divergence 
encountered in i3TI . 



5.2.2 Divergences in the Total Energy 

In this subsection we are going to use heat-kernel knowledge to determine the di- 
vergence structure in the total energy. We consider a general cylinder of the type 
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'^ — RxY, where Y is an arbitrary smooth two dimensional region rather than 
merely being the disc. As a metric we have ds^ = dz^ + dY^ from which we obtain 
that the zeta function (density) associated with the Laplacian on ^ is {Res > 3/2) 

Ay ^ ' Ay 



CyIs--]. (163) 



2k r{s) "' y 2 

Here Xy are the eigenvalues of the Laplacian on Y, and ^y (■*) is the zeta function 
associated with these eigenvalues. In the zeta-function scheme the Casimir energy 
is defined as 



£cas - ^M^' C ^ 



(164) 

s=0 



which, in the present setting, turns into 



Expanding this expression about s = Q, one obtains 



(165) 



£cas = ^Cy(-l) + 5^ (CF(-l)[21n(2M)-l] + C;(-l)) + €?(.). (166) 

The contribution associated with tl,y{—\) can be determined solely from the heat- 
kernel coefficient knowledge, namely 

Cf(-1) = -«4, (167) 

in terms of the standard 4th heat-kernel coefficient. The contribution coming from 
C,y{—\) can in general not be determined. But as we see, at least the divergent term 
can be determined entirely by the heat-kernel coefficient. 

The situation considered in the Casimir energy calculation is a 5-function shell 
along some smooth line I, in the plane (here, a circle of radius a). The manifolds 
considered are the cylinder created by the region inside of the line, and the region 
outside of the line; from the results the contribution from free Minkowski space has 
to be subtracted to avoid trivial volume divergences (the representation in terms of 
the Bessel functions already has Minkowski space contributions subtracted). The 5- 
function shell generates a jump in the normal derivative of the eigenf unctions; call 
the jump U (here, U = X/a). The leading heat-kernel coefficients for this situation, 
namely for functions which are continuous across the boundary but which have a 
jump of the first normal derivative at the boundary, have been determined in 1891 ; 
the relevant 04 coefficient is given in Theorem 7.1, p. 139 of that reference. The 
results there are very general; for our purpose there is exactly one term that survives, 
namely 
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(168) 



a4 = 


24a: Je 


'div 


^ f iMT''^ 


Cas 


\92nhJE ^ 



which shows that 

(169) 

So no matter along which Hne the 5-function shell is concentrated, the first two 
orders in a weak-coupling expansion do not contribute any divergences in the total 
energy. But the third order does, and the divergence is given above. 

For the example considered, as mentioned, U — X/a is constant, and the inte- 
gration leads to the length of the line which is 2na. Thus we get for this particular 
example 

1 A3 



'^cZ = ^^^-j- (170) 



^div 

967ZS fl2 

[Compare this with the corresponding divergence for a sphere, dSSTl.l This can be 
easily checked from the explicit representation we have for the energy. We have 
already seen that the first two orders in X identically vanish, while the part of the 
third order that potentially contributes a divergent piece is 

The m = contribution is well behaved about s = Q; while for the remaining sum 
using 

Kl{mz)li{mz) ^ ^,^^^^^, (172) 

we see that the leading contribution is 

^(3) ^ f ^2-2.s r^,,2-2.v d 1 1 



127rfl2^,^j JQ ~~ dz8m3 (1+72)3/2 

l^ „ ,. . . r . ,_,.. d 1 



_2-2s 



^«.+2,n^l)£ii±i).^,.,,.), ,,,3, 



96na^^^ ' T{l,ll) 96nah 

in perfect agreement with the heat-kernel prediction ( I170l i. 

5.3 Strong Coupling 

The strong-coupling limit of the energy ( 1149b . that is, the Casimir energy of a Dirich- 
let cylinder, 

^^ = -J-^ E / dxx^-^\nUx)K,„{x), (174) 

ona^ ,„r^^jQ ax 
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was worked out to high accuracy by Gosdzinsky and Romeo 

^r, 0.000614794033 

© = ;; . 



(175) 



It was later redone with less accuracy by Nesterenko and Pirozhenko 119011 . 

For completeness, let us sketch the evaluation here. We carry out a numerical cal- 
culation (very similar to that of 1901 ) in the spirit of Sect. l5.2.TI We add and subtract 
the leading uniform asymptotic expansion (for m = Q the asymptotic behavior) as 
follows: 



^« 



1 



8;rfl2 



— 2 / dxx 



ln(2xIo{x)Koix))- 



1 1 



l+x2 



m=l 



2V / dxx^- 
hh Ax 



In 



{2xI,n{x)K,n{x))~\n(-)~\'^ 



2|-+V I / djfjc^ — ln2x + 2 y / dxx^— Inxf 
U h Jo dx ,tJo dx 




(176) 



In the first two terms we have subtracted the leading asymptotic behavior so the 
resulting integrals are convergent. Those terms are restored in the fourth, fifth, and 
sixth terms. The most divergent part of the Bessel functions are removed by the in- 
sertion of 2x in the corresponding integral, and its removal in the third term. (As 
we've seen above, such terms have been referred to as "contact terms," because if a 
time-spUtting regulator, e'^^, is inserted into the frequency integral, a term propor- 
tional to 5(t) appears, which is zero as long as T 7^ 0.) The terms involving Bessel 
functions are evaluated numerically, where it is observed that the asymptotic value 
of the summand (for large m) in the second term is I /32m". The fourth term is 
evaluated by writing it as 



21im y 



„2-.s 



s^O 



,1-s 



l+t 



= 2C'(-2) = - 



2;r2' 



(177) 



while the same argument, as anticipated, shows that the third "contact" term is zerojj 
while the sixth term is 



^ This argument is a bit suspect, since the analytic continuation that defines the integrals has no 
common region of existence. Thus the arguinent in the following subsection may be preferable. 
However, since that terin is properly a contact term, it should in any event be spurious. 
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1 
lim 

2 .v^o 



'as)^\ 



- = -\n2n. (178) 

s 4 



The fifth term is elementary. The result then is 



^«= ^ 



- — ^ (0.010963 - 0.0227032 + + 0.0304485 + 0.21875 - 0.229735) 

1 , , 0.0006146 

(0.007724) = , (179) 



which agrees with (1175b to the fourth significant figure. 

5.3.1 Exponential Regulator 

As in Sect. 15.2.11 it may seem more satisfactory to insert an exponential regulator 
rather than use analytic regularization. Now it is the third, fourth, and sixth terms in 
( |176t that must be treated. The latter is just the negative of ( |161t . We can combine 
the third and fourth terms to give using ( II6OI 1 

' 2 ^^ '^-^ '^^ ' (180) 



52 5270 7j-A,?P-A{?-e-V 

The latter integral may be evaluated by writing it as an integral along the entire z 
axis, and closing the contour in the upper half plane, thereby encircling the poles at 
i5 and at 2inn, where n is a positive integer. The residue theorem then gives for that 
integral 

2^^ ^(^) ri8n 

"5?"2^' ^^^^^ 

so once again we obtain the same finite part as in ( I177l i. In this way of proceeding, 
then, in addition to the finite part in ( |179t , we obtain divergent terms 

'^'''^"64^^8^^^252+^^253' (182) 

which, with the previous interpretation for 5, implies divergent terms in the energy 
proportional to Lja (shape), L (length), and ah (area), respectively. Such terms pre- 
sumably are to be subsumed in a renormalization of parameters in the model. Had a 
logarithmic divergence occurred [as does occur in weak coupling in 6()\?y\ such a 
renormalization would apparently be impossible — however, see ll37l . 



5.4 Local Energy Density 

We compute the energy density from the stress tensor (|145l l, or 
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(r00) = ±(50a0' + V.V')G(x,x') 



1 



l6nH J- 



dk 



dco ^ 



■^V^G(x,x) 



(o'^ + k^ + ^- + drd,, ] g{r, r') 



2^-drrdrg{r,r) 



(183) 



We omit the free part of the Green's function, since that corresponds to the energy 
that would be present in the vacuum in the absence of the cyHnder When we insert 
the remainder of the Green's function ( I144bb . we obtain the following expression 
for the energy density outside the cylindrical shell: 



„(r) = (ro«-rp 



2co^+K- 



r > a. 



A 



I6n^ J- 



'^L'\?^-i 



ll{Ka) 



+ Xl,„{Ka)K,„{Ka) 
Ki{Kr) + K^K'^^[Kr)-2^~r^Kl{Kr) 



(184) 



The factor in square brackets can be easily seen to be, from the modified Bessel 
equation, 



, ,, , 1-4^ Id d ., , 



„,„,,..,. (185) 

2 r or dr 

For the interior region, r <a,v4C have the corresponding expression for the energy 
density with /,„ o K,„. 



5.5 Total and Surface Energy 

We first need to verify that we recover the expression for the energy found in 
Sect. 15.11 So let us integrate expression ( |184t over the region exterior of the cylinder, 
and the corresponding interior expression over the inside region. The second term 
in (1185b is a total derivative, while the first is exactly the one evaluated in Sec. 15. II 
The result is 



2n / drru{r) 
Jo 



1 y 


I 


(l-4^)l 


A 



dxx^—ln[l+Unix)K„{x)] 



4na^ Jo 



djc 



dxx 



Jt 



I,n{x)Km{x) 
+ Xl,n{x)Kn{x)' 



(186) 



The first term is the total energy ( 11491 ), but what do we make of the second term? In 
strong coupUng, it would represent a constant that should have no physical signifi- 
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cance (a contact term — it is independent of a if we revert to the physical variable k 
as the integration variable). In general, however, there is another contribution to the 
total energy, residing precisely on the singular surface. This surface energy is given 
in general by IgOl SB 121 23 



1 - 4£ r 
E = ^j,dS-VG{x,x'] 



(187) 



as given for ^ = in ( |20a| ), where the normal to the surface is out of the region in 
question. In this case it is easy to see that E exactly equals the negative of the second 
term in (II86I 1. This is an example of the general theorem (I2TI) 



(dr)M(r) +£■ = £, 



(188) 



that is, the total energy E is the sum of the integrated local energy density and the 
surface energy. The generalization of this theorem, ( |187t and ( |188t , to curved space 
is given in ||57l . A consequence of this theorem is that the total energy, unlike the 
local energy density, is independent of the conformal parameter ^ . (Note that this 
surface energy vanishes when ^ = 1/4 as Fulling has stressed |[93l .) 



5.6 Surface Divergences 



We now turn to an examination of the behavior of the local energy density (11841 ) as 
r approaches a from outside the cylinder To do this we use the uniform asymptotic 
expansion dSTT i. Let us begin by considering the strong-coupling limit, a Dirichlet 
cyhnder. If we stop with only the leading asymptotic behavior, we obtain the ex- 
pression 



u{r) 



1 



in= Jo 



— J / dKKi y 



.-n'X , 



-{1-^)k' 



n 1 

2mt 7} 



-)c" + (l-4<^)hf" + ^ 



(A^oc), 



nt 
2m 



(189) 



where 



X 



T7(z)-TJ Z 



(190) 



and we have replaced the integral over k and ^ by one over the polar variable K as 
before. Here we ignore the difference between r and a except in the exponent, and 
we now replace K by mz/a. Close to the surface. 



2 r — a 

X~- , r-a-^r, 

t r 



(191) 
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and we carry out the sum over m according to 










2f; mV"'^~ 2// - 


12 


3 


r n\A- 


(192) 



Then the energy density behaves, as r ^ «+, 

3 1 



^ ^ r(l-6^). (193) 



16%"^ (r— I 

This is the universal surface divergence first discovered by Deutsch and Candelas 
lIMI and seen for the sphere in ( I127ab i74|. It therefore occurs, with precisely the 
same numerical coefficient, near a Dirichlet plate ||36l . Unless gravity is considered, 
it is utterly without physical significance, and may be eliminated with the conformal 
choice for the parameter (^ , ^ = 1 /6. 

We will henceforth make this conformal choice. Then the leading divergence 
depends upon the curvature. This was also worked out by Deutsch and Candelas 
||6T| ; for the case of a cylinder, that result is 

/20n^ r[r — ay 

exactly 1/2 that for a Dirichlet sphere of radius a (I130I I ITmI , as anticipated from 
the general analysis summarized in ( |59] l. Here, this result may be straightforwardly 
derived by keeping the 1 /m corrections in the uniform asymptotic expansion (ISTT l. 
as well as the next term in the expansion of x^ (11281 ). 



5.6.1 Weak Coupling 

Let us now expand the energy density ( |184| l for small coupling, 

"W = -T^ r ^^ r '^^ t ll{Ka)l^{-Xril{Kc,)Kl{Ka) 

iOTT^ 7-00 7-00 „,=_„ „^o 



->c2 + (l-4^)(^2 + ^ 



(195) 

If we again use the leading uniform asymptotic expansions for the Bessel functions, 
we obtain the expression for the leading behavior of the term of order X", 

"'"^W - ^ f-^y C^- i m3-"e-"'^r«-'(f2 + l -8^). (196) 

m=l 



87r2r4 V 2/70 
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The sum on m is asymptotic to 

CO / f \ ^^" 

so the most singular behavior of the order X" term is, as r ^- a+. 

This is exactly the result found for the weak-coupling limit for a 5-sphere ( |127b| l 
ITMII and for a 5-plane ( |48l ) P31 . so this is also a universal result, without physical 
significance. It may be made to vanish by choosing the conformal value ^ = 1/6. 

With this conformal choice, once again we must expand to higher order. We use 
the corrections noted above, in dSTl i and (1128b . (|129t . Then again a quite simple 
calculation gives 

«w~f 3N«l!i_l)(!i±^)n3_«) ,,QQ. 

" ^^^> 28807rV+i('-«)^-"' ' ^ ^ 

which is analytically continued from the region 1 < Ren < 3. Remarkably, this 
is exactly one-half the result found in the same weak-coupling expansion for the 
leading conformal divergence outside a sphere dlBll l fMl . Therefore, like the strong- 
coupling result (|194| i, this limit is universal, depending on the sum of the principal 
curvatures of the interface. 

In ||37l we considered a annular shell of finite thickness, which as the thickness 
5 tended to zero gave a finite residual energy in the annulus, in terms of the energy 
density u in the annulus, 

^ann-2;r5a«^(l-4^)-A^ f H dKaKa-^^^^^^^^^^^^^ ^i, (200) 
4na^ ,„t^„jQ I + Xl,j,{Ka)K,„{Ka) 

which is exactly the form of the surface energy given by the negative of the second 
term in (1186b . In particular, note that the term in S' of order X^ is, for the conformal 
value ^ = 1/6, exactly equal to that term in the total energy S' ( 11491 ): [see ( 11711 )1 

#(3'=(f(^l (201) 

This means that the divergence encountered in the global energy ( 1170b is exactly 
accounted for by the divergence in the surface energy, which would seem to provide 
strong evidence in favor of the renormalizablity of that divergence. 
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6 Gravitational acceleration of Casimir energy 

We will here show that a body undergoing uniform acceleration (hyperbolic motion) 
imparts the same acceleration to the quantum vacuum energy associated with this 
body. This is consistent with the equivalence principle that states that all forms of 
energy should gravitate equally. A general variational argument, which, however, 
did not deal with the divergent parts of the energy, was given in 1221 . This section is 
based on 



6.1 Green's Functions in Rindler Coordinates 

Relativistically, uniform acceleration is described by hyperbolic motion, 

z^^coshT and f = ^sinhT. (202) 

Here the proper acceleration of the particle described by these equations is (^ ^ ' , and 
we have chosen coordinates so that at time f = 0, z(0) = ^. Here we are going to 
consider the corresponding metric 

As^ == -df 2 + dz^ + Ax- + Ay- = -^ -Ai^ + A^- + Ax^ + d/. (203) 

In these coordinates, the d' Alembertian operator takes on cylindrical form 



where _L refers to the x-y plane. 



6.1.1 Green's Function for One Plate 

For a scalar field in these coordinates, subject to a potential V (x), the action is 

W ^ JA\^-g{x)^{^{x)), (205) 

where x = (T,x,y, i^) represents the coordinates, d'*x = drd^djcdy is the coordi- 
nate volume element, g^iv{x) = diag(—^^, + 1,-1-1, +1) defines the metric, g{x) = 
Ae.igfxv (x) = — ^ ' is the determinant of the metric, and the Lagrangian density is 

■^(•^ W) = -\gi,v{x)d^<^{x)d',^{x) - \y{x)<^{xf, (206) 

where for a single semitransparent plate located at B,i 
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y(x) = A5(^-^i), 
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(207) 



and A > is the coupling constant having dimensions of mass. More explicitly we 
have 



W = (fx^ 



my-m'-^^^^>'--'>^' 



(208) 



Stationarity of the action under an arbitrary variation in the field leads to the equa- 
tion of motion 



1 ^2 I d ^ d 



-V{x) 



(j){x)=0. 



The corresponding Green's function satisfies the differential equation 



(209) 



1 d^ 1 <9 , <9 . 



G{x,x') = ^^^y^5(T-T')5(xx -x'^). 

(210) 

Since in our case V{x) has only £, dependence we can write this in terms of the 
reduced Green's function g{£, ,£,'), 



where §(i^,i^') satisfies 









(^,r)^^(^-^') 



^ 



(212) 



We recognize this equation as defining the semitransparent cylinder problem dis- 
cussed in Sect. |5] 1371 . with the replacements 



fl— >^l, m — > ^ = — /o, iC— >fc = A:x, A — > A(^i 



(213) 



so that from (I144ab and (I144bl l we may immediately write down the solution in 
terms of modified Bessel functions. 



^(<^,^')-/f(^^<)^c(^^> 



i+A^i/^(^^i)/:f(fe^i) ' 



^,r<<^i 



(214a) 






(214b) 



Note that in the strong-coupling limit, A — > oo, this reduces to the Green's function 
satisfying Dirichlet boundary conditions at i^ = ^i ■ 
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6.1.2 Minkowski-space Limit 

To recover the Minkowski-space Green's function for the semitransparent plate, we 
use the uniform asymptotic expansion (Debye expansion), based on the Hmit 

i^^oo, ^i^°o, ^-^1 finite, C^°°i C/^i finite . (215) 

For large C, we use dSTT i with x = ^z = A:^ , for example. Expanding the above expres- 
sions (|214a| i, ( |214bt around some arbitrary point ^o. chosen such that the differences 
^ ~ £,0, £,' — £,(), and £,i — £,o are finite, we find for the leading term, for example, 

,/Wli;mKi;ik^') - :^e^(«-«'), (216) 

where K^ = k^ + t,^, t, = ^/^q. In this way, taking for simplicity t,Q = t,i, we find 
the Green's function for a single plate in Minkowski space. 



6.1.3 Green's Function for Two Parallel Plates 

For two semitransparent plates perpendicular to the t^-axis and located at i^i, B,2, 
with couplings X\ and X2, respectively, we find the following form for the Green's 
function: 

5(Si? J - ^<"^> ^ "': 

^,r<^i, (218a) 

^ ^ Xx^,l\ + X2^2ll + hhili2hh{hKl-hK2)^^ 
— i<J^> T J^J^h 

^,r>fe, (218b) 

X2^2Klil+h^iKih) 

A 

A A 

^i<^,r<fe, (218c) 



where 



A = {l+Xi^lKih){l+X2^2K2l2)-^l^2^l^2ltK^, (219) 

and we have used the abbreviations /i = Ir{k^i), I = Ir{k^), I, = Ir{k^'), etc. 

Again we can check that these formulas reduce to the well-known Minkowski- 
space limits. In the ^0 ^ °° limit, the uniform asymptotic expansion (ISTl i gives, for 

^l<^,^'<^2 
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A1A2 , 



^ogi^A') ^ g'^'K^,^') = ^e-^i^-^'i + ^ 



4K-2 



2coshK-(<^-^') 



2k- V 2k-/ 2k- V 2k-/ 

(220) 

where (a = ^2 - ^1) 

which is exactly the expected resuh ( fTal i. dHJ. The correct Hmit is also obtained in 
the other two regions. 



6.2 Gravitational Acceleration of Casimir Apparatus 

We next consider the situation when the plates are forced to "move rigidly" ||94l in 
such a way that the proper distance between the plates is preserved. This is achieved 
if the two plates move with different but constant proper accelerations. 

The canonical energy-momentum or stress tensor derived from the action ( 1205b 
is 

Ta^{x) = daii>{x)dp(i>{x)+gafi{x)^{(i>{x)), (222) 

where the Lagrange density includes the 5-function potential. The components re- 
ferring to the pressure and the energy density are 



-\m)'^\m-\^^^'^'-T 



J-MW = iT^(S) +ilS) --A^^^f-^y{')^\ (223a) 



2 



1 ^ ..A _^ ^ (d^\\i fd^y , 1 ,„ ^,2 , i,,,..x.,2 



The latter may be written in an alternative convenient form using the equations of 
motion ( |209t : 

1 fd^y^ 1 '5' . id / , d \ ^2 



which is the appropriate version of ( fT9l ) here. The force density is given by 

h = -^dy{^^T\) + \T^'dxgf,v, (225) 

\ -g 2 



or in Rindler coordinates 
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/^=-l^^(^r«^)-^r«o. (226) 

When we integrate over all space to get the force, the first term is a surface term 
which does not contributeO 

^ = |d^^/^ = -|||roo. (227) 

This could be termed the Rindler coordinate force per area, defined as the change 
in momentum per unit Rindler coordinate time T per unit cross-sectional area. If 
we multiply ^ by the gravitational acceleration g we obtain the gravitational force 
per area on the Casimir energy. This result ( I227l i seems entirely consistent with the 
equivalence principle, since ^ "^Too is the energy density. Using the expression (I2241 i 
for the energy density, taking the vacuum expectation value, and rescaling i^ = i^^, 
we see that the gravitational force per cross sectional area is merely 

^ = |d^^|^CV^,^). (228) 

This result for the energy contained in the force equation ( 1228b is an immediate 
consequence of the general formula for the Casimir energy ( fT6] l |[38J. 

Alternatively, we can start from the following formula for the force density for a 
single semitransparent plate, following directly from the equations of motion ( 12091 ), 

/^ = i(/.25^A5(^-^i). (229) 

The vacuum expectation value of this yields the force in terms of the Green's func- 
tion. 



(230) 



6.2.1 Gravitational Force on a Single Plate 

For example, the force on a single plate at ^\ is given by 



* Note that in previous works, such as I45II461 . the surface term was included, because the integra- 
tion was carried out only over the interior and exterior regions. Here we integrate over the surface 
as well, so the additional so-called surface energy is automatically included. This is described in 
the argument leading to j20at . Note, however, if i226t is integrated over a small interval enclosing 
the 5-function potential, 

where AT^^ is the discontinuity in the normal-normal component of the stress density. Dividing 
this expression by ^\ gives the usual expression for the force on the plate. 
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1 /-dCd^k, 

2 J llnf' 



^-'\^ I 7^1n[l+^^i^?(^^i)^?(^^i)]' (231) 



Expanding this about some arbitrary point £,q, with ^ — ^£,0, using the uniform 
asymptotic expansion (ISTT i. we get {k^ = k^ + t,"^) 

From this, if we introduce polar coordinates for the k- 1^ integration, the coordinate 
force is 



(233) 



1 r dyy^ 



96n^a^ Jo 1+y/Xa 
where for example 



„2 a ..2 1 



(C) = -/ dcos0cos^0K-' = -K-^ (234) 



The divergent expression (I233l l is just the negative of the quantum vacuum energy 
of a single plate, seen in ( fTTl i and (l43l l. 



6.2.2 Parallel Plates Falling in a Constant Gravitational Field 

In general, we have two alternative forms for the gravitational force on the two-plate 
system: 

^ = -(^^,+%)^/f|?ln4- (235) 

A given in ( I219l l. which is equivalent to (1228b . (In the latter, however, bulk energy, 
present if no plates are present, must be omitted.) From either of the above two 
methods, we find the coordinate force is given by 

^ = --— 5-/ dfCK-'lnzio, (236) 

4?:^^ Jo 

where Aq = e^^'^"A, A given in (1221b . The integral may be easily shown to be 

, ,„ 1 J 1 J 1 

i / . 3 y+^i y+^i 



96nVJo ^'' (^ + i)(^ + il,v_i 
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1 



96nV 







dy/ 



aTci aqTi 



(237a) 



= -(^, + 4/1+4/2), (237b) 

which is just the negative of the Casimir energy of the two semitransparent plates 
including the divergent pieces — See ( fTTI i 1451 |46l . Note that (f^/,, / = 1,2, are simply 
the divergent energies ( |233l l associated with a single plate. 



6.2.3 Renormalization 

The divergent terms in ( |237b| l simply renormalize the masses (per unit area) of each 
plate: 

-Etotai = mi+m2 + S',n + ^42 + Sc 

= Mx+M2 + Sc, (238) 

where m,- is the bare mass of each plate, and the renormalized mass is M, = m,- + (f^/,. 
Thus the gravitational force on the entire apparatus obeys the equivalence principle 

g^ = -g{Mi +M2 + Sc)- (239) 

The minus sign reflects the downward acceleration of gravity on the surface of the 
earth. Note here that the Casimir interaction energy S^ is negative, so it reduces the 
gravitational attraction of the system. 



6.3 Summary 

We have found, in conformation with the result given in 122], an extremely simple 
answer to the question of how Casimir energy accelerates in a weak gravitational 
field: Just like any other form of energy, the gravitational force F divided by the 
area of the plates is 

J = -gSc (240) 

This is the result expected by the equivalence principle, but is in contradiction 
to some earlier disparate claims in the literature ||95l l96l l97l l98l l99l . Bimonte et 
al. II 1001 now agree completely with our conclusions. This result perfectly agrees 
with that found by Saharian et al. IIOII for Dirichlet, Neumann, and perfectly con- 
ducting plates for the finite Casimir interaction energy. The acceleration of Dirichlet 
plates follows from our result when the strong coupling limit A ^ 00 is taken. What 
makes our conclusion particularly interesting is that it refers not only to the finite 
part of the Casimir interaction energy between semitransparent plates, but to the di- 
vergent parts as well, which are seen to simply renormalize the gravitational mass 
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of each plate, as they would the inertial mass. The reader may object that by equat- 
ing gravitational force with uniform acceleration we have built in the equivalence 
principle, and so does any procedure based on Einstein's equations; but the real non- 
triviality here is that quantum fluctuations obey the same universal law. The reader 
is also referred to the important work on this subject by Jaekel and Reynaud 01021 . 
and extensive references therein. 



7 Conclusions 

In this review, I have illustrated the issues involved in calculating self-energies in the 
simple context of massless scalar fields interacting with 5-function potentials, so- 
called semi transparent boundaries. This is not as unrealistic as it might sound, since 
in the strong coupling limit this yields Dirichlet boundary conditions, and by using 
derivative of 5-function boundaries, we can recover Neumann boundary conditions. 
Thus, where the boundaries admit the separation into TE and TM modes, we can 
recover perfect-conductor boundaries imposed on electromagnetic fields. 

We have examined both divergences occurring in the total energy, and diver- 
gences which appear in the local energy density as boundaries are approached. The 
latter divergences often have little to do with the former, because the local diver- 
gences may cancel across the boundaries, and they typically depend on the form 
(canonical or conformal, for example) of the local stress-energy tensor The global 
divergences apparently can always be uniquely isolated, leaving a unique finite self- 
energy; in some cases at least the divergent parts can be absorbed into a renormal- 
ization of properties of the boundaries, such as their masses. It is expected that if the 
ideal boundaries were represented as a solitonic structure arising from a background 
field, this "renormalization" idea could be put on a more rigorous footing. 

Evidence for the consistency of this view occurs in the parallel plate configura- 
tion, where we show that the finite interaction energy and the divergent self-energies 
of each plate exhibit the same inertial and gravitational properties, that is, are each 
consistent with the equivalence principle. Thus it is indeed consistent to absorb the 
self-energies into the masses of each plate. We hope to prove in the future that this 
renormalization consistency is a general feature. 

In spite of the length of this review, we have barely scratched the surface. In 
particular, we have not discussed how the divergent contributions of the local stress 
tensor are consistent with Einstein's equations II103I . We have also only discussed 
simple separable geometries, where the equations for the Green's functions can be 
solved on both the inside and the outside of the boundaries. This excludes the exten- 
sive work on rectangular cavities, where only the sum over interior eigenvalues can 
be carried out lfT6l[T7l[T8l[T9l[T04l . There are some numerical coincidences, for ex- 
ample between the energy for a sphere and a cube, but since divergences have been 
simply omitted by zeta-function regularization, the significance of the latter results 
remains unclear There are a few other examples where the interior Casimir con- 
tribution can be computed exactly, while the exterior problem cannot be solved, an 
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example being a cylinder with cross section of an equilateral triangle. Such results 
seem more problematic than those we have discussed here. 

We also have not discussed semiclassical and numerical techniques. For exam- 
ple, there is the extremely interesting work of Schaden II105L who computes a very 
accurate approximation for the Casimir energy of a spherical shell using optical path 
techniques. The same technique gives zero for the cylindrical shell, not the attractive 
value found in |27;|, which is not surprising. Not unrelated to this technique is the 
exact worldline method of Gies and collaborators II106II1071I108II . which is able to 
capture edge effects. The optical path work of Scardicchio and Jaffe ||109llll0llllll 
should also be cited, although it is largely restricted to examining the forces between 
distinct bodies. This review also does not refer to the remarkable progress in numer- 
ical techniques, some of which are related to the multiple scattering approach — for 
some recent references see II112II113JI . and the contributions to this volume by Emig, 
Jaffe, and Rahi and by Johnson — , which however, have not yet been turned to ex- 
amining self-interactions. 

The central issue is the meaning of Casimir self-energy, and how, in principle, it 
might be observed. Probably the right direction to address such issues is in terms of 
quantum corrections to solitons — for example, see 11 14111 15lll 1611 . The issues being 
considered go to the very heart of renormalized quantum field theory, and likely to 
the meaning and origin of mass, a subject about which we in fact know very little. 
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